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Majorana representation of quantum states by a constellation of n 'stars' (points on the sphere) 
can be used to describe any pure state of a simple system of dimension n + 1 or a permutation 
symmetric pure state of a composite system consisting of n qubits. We analyze the variance of the 
distribution of the stars, which can serve as a measure of the degree of non-coherence for simple 
systems, or an entanglement measure for composed systems. Dynamics of the Majorana points 
induced by a unitary dynamics of the pure state is investigated. 
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I. INTRODUCTION 

Non-classical correlations between composite quantum 
systems became a subject of an intense current research. 
A particular kind of such correlations, called quantum en- 
tanglement attracts special attention of theoretical and 
experimental physicists - see [l| and references therein. 
Entanglement in quantum systems consisting of two sub- 
systems is nowdays relatively well understood, but sev- 
eral significant questions concerning multiparticle sys- 
tems remain still open. 

One of the key issues is to workout a practical entangle- 
ment measure, used to quantify the quantum resources 
of a given state. Although various measures of quantum 
entanglement are known [H-[3|; they are usually difficult 
to compute. An important class of entanglement mea- 
sures can be formulated within the geometric approach 
to the problem 0, Q . For a given quantum state one can 
study its minimal distance to the set of separable states. 
Various distances Q can be used for this purpose and 
the minimization can be performed with respect to the 
mixed [f| or pure @, H| separable states. 

Working with pure quantum states of a n-qubit system 
it is possible to distinguish a class of states symmetric 
with respect to permutations of all subsystems. This 
class of symmetric pure quantum states can be identified 
with the set of all states of a single system described in 
a N = n + 1 dimensional Hilbert space. 

Analyzing the space of pure states belonging to N di- 
mensional Hilbert space it is useful to distinguish the 
class of spin coherent states. These states \8, 0), labeled 
by a point on the sphere can be obtained by action of 
the Wigner rotation matrix R$,$ on the maximal weight 
state \j,j), so they are also called SU(2) coherent states. 
Here j = n/2 is the maximal eigenvalue of the compo- 
nent J z of the angular momentum operator. Making use 
of the stereographic projection one can map the sphere 
into the plane. The coherent states are then labeled by a 
complex number a and their expansion in the eigenbasis 
of J z reads UGH 



\ip) = X)fc=o Cfc l^)- -^ can b e expanded in the coherent 
states representation, 
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Consider now an arbitrary state of an TV-level system, 



The function Q^{a) is called the Husimi function (or Q— 
function) of the state 1^) and it can be interpreted as a 
probability density on the plane. As it can be associated 
with a polynomial of order n of a complex argument, 
it is uniquely represented by the set if its n roots Zi, 
i = 1, ... n, which may be degenerated. With help of the 
inverse stereographic projection one can map these points 
back on the sphere. The collection of these n points on 
the sphere, corresponding to the zeros of the Husimi func- 
tion ([2]), represents uniquely the state € rl n +i- This 
approach of Majorana and Penrose [12| leads to the 
so-called stellar representation of a state, as each Ma- 
jorana point on the sphere corresponding to 0j, can be 
interpreted as a star on the sky. 

For any coherent state |cv) = \9, <fi) all n stars sit in 
a single point antipodal to the vector pointing in the 
direction (8, <fi). If the Majorana points are located in a 
neighborhood of a single point, the corresponding state 
is close to be coherent. In contrast, a generic random 
state, for which the distribution of stars is uniform on 
the sphere [HI, [3 is far from being coherent. The degree 
of non-coherence of \ip) can be characterized e.g. by its 
Fubini-Study distance to the closest coherent state, or by 
the Monge distance equal to the total geodesic distance 
on the sphere all stars have to travel to meet in a single 
point jl 51 ] . 

The stellar representation, originally used for states of 
a simple system of n + 1 levels, can be also used to ana- 
lyze composite systems consisting of n qubits, under the 
assumption that the states are symmetric with respect 
to permutations of all subsystems f2ljj . Due to this 
symmetry the investigation of such states becomes easier 
and estimation of their geometric measure of entangle- 
ment can be simplified [22|. 

Note that any separable state of the n-qubit system 
can be considered as a coherent state with respect to the 
composed group [SU(2)] n . Thus entangled states of a 
composite system correspond to non-coherent states of 
the simple system with n + 1 levels [23| , while the degree 
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of entanglement can be identified with the degree of non- 
coherence. 

In this work we propose to characterize pure quantum 
states by the position of the barycenter of its Majorana 
representation. We will then describe non-classical prop- 
erties of a state by the average and the variance of the 
corresponding distribution of the Majorana points. The 
latter quantity has a simple geometric interpretation as 
a function of the radius of the barycenter of the stars, 
which is situated inside the Bloch ball. The same quan- 
tity can thus be applied to characterize the degree of 
non-coherence of a pure state of an (n + l)-level sys- 
tem, or simultaneously, as a degree of entanglement for 
the corresponding permutation symmetric pure state of 
n-qubit system. 

In a similar way the measure of quantumness of a state 
of n + 1 dimensional system introduced by Giraud et 
al. [24| can be used to characterize the entanglement of 
permutationally symmetric pure states of n qubits. We 
are also going to study time evolution of quantum states 
and an associated dynamics of stars. In particular we 
investigate how the barycentric measure of entanglement 
varies with time. 

This work is organized as follows. In section II we 
review the Majorana representation used for symmetric 
pure states of a multiqubit quantum system. In Section 
III the position of the barycenter of Majorana points and 
their variance is investigated, as it may serve as a charac- 
terization of quantum entanglement. A family of n-qubit 
states maximally entangled with respect to this measure 
is studied in section IV. Unitary dynamics of quantum 
states and the associated dynamics of Majorana points 
on the sphere is analyzed in Section V. 



II. MAJORANA REPRESENTATION OF 
MULTI-QUBIT STATES 

Consider a n-qubit pure state symmetric with re- 
spect to permutation of subsystems: 

I ^ = 77f E I ^(l)> I &r(2)) ' ' ' I &r(70>- (3) 

The sum is taken over all permutations ir, while the nor- 
malization constant reads 

n 

K = nl EIIi^i^))!- 

Any one-qubit pure state can be represented as | <f>i) — 
cos f | 0) + e^' sin f | 1). Making us of the the Dicke 
states | S ni k) with k = 0, 1, . . . ,n [25j |. 

_ i 

\ s ^ = (t) 2 T, p p\°^ n ' k \ 1 ^^ ( 4 ) 

^ ' p 



we can represent the state \ip) as their superposition, 

i^> = X>is n , fc > = x>( n k ) 'j2p P \or n - k \^ k - 

k=0 k=0 ^ ' p 

(5) 

Here dk denote complex coefficients, (^) is the binomial 
coefficient and the sum ^ P p goes over all states of n- 

p 

qubits with exactly k qubits in the state |1) and n — k 
qubits in the state |0). It comes out that quotients of the 

coefficients .™ s ." T t) . = e~ l$i cot(6>i/2) determining the 

e i sin 

one-qubit state \4>i) can be obtained as roots Zi of the 
polynomial 

p(z) = ]>>(-i) fe r) z k . (6) 

The polynomial defined above has degree D < n and 
consequently as many roots. In case of D < n the remain- 
ing (n — D) coefficients cos If we set to unity [l6j . 

In this way we can represent any n-qubit state sym- 
metric with respect to permutations as n points on the 
Bloch sphere related to spin-i (1-qubit) states | <j)i) = 
cos f- | 0) + e 4$i sin y | I). These n points on the Bloch 
sphere define a stellar representation of a state and are 
called stars or Majorana points (MP), while polynomial 
((B]) is called the Majorana polynomial [Tll |. 

Note that the same constellation of n stars determines 
on one hand a pure state of the simple system described 
in TV = n + 1 dimensional Hilbert space [USUI- On the 
other hand, it describes also a permutationally symmetric 
pure state of a n qubit system [l?], H3, [H| and belongs 
to the Hilbert space of dimension 2". 

To establish a one-to-one link between both problems 
it is sufficient to identify bases in both spaces. In the 
case of a simple system described in the Hilbert space of 
size n+ 1 we select the standard eigenbasis of the angular 
momentum operator J z , 

\j,m), m = -j,...,j (7) 

with j = n/2. A state from this basis is represented by 
a constellation of n/2 + m stars at the north pole and 
the remaining n/2 — m stars at the south pole as shown 
in Fig. [TJ Thus the state | j, m) € Hn can be identified 
with the Dicke state \S n , k ) € Uf n defined in flU) with 
k = n/2 — to. In this way any constellation of stars can 
be used to describe pure quantum state of two different 
physical systems: an (n+l)-level system and a symmetric 
state of n qubits. In following sections we show that this 
link can be extended also for quantum dynamics. Any 
discrete unitary dynamics of an n-qubit system, which 
preserves the permutation symmetry, can be represented 
by an unitary matrix of size 2™ with a block structure in 
the computational basis. The diagonal block of size n + 1 
is unitary and it defines the corresponding dynamics of 
a simple quantum system with n + I levels. 
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simple composite 
systems systems, n=2 

|j,m> 



lu> 



1,0> 



1,-1> 



|00> 

[|10>+|01>] 

lll> 



simple composite 
systems systems, n=3 



3/2,3/2> 



000 > 



3/2,l/2> ( \ [|001>+|010>+[100>] 



|3/2,-l/2>^_J [|U0>+|101>+|011>] 
T 

|3/2,-3/2>^^ | m> 



10) 




FIG. 1. Stellar representation of the orthogonal basis in (n+1) 
dimensional Hilbert space describing a simple system (states 
\j, m) with j — n/2) with the corresponding basis in the sub- 
space of permutation symmetric states of n-qubit systems 
plotted for n — 2 and n = 3. 



FIG. 2. (color online) Majorana representation of a symmetric 
2-qubit state of the form = |0) and \<f>2) = cos||0) + 
sin with 8 = ~7r. Majorana points are shown as red and 
green balls. The black cube denotes the barycenter of MP 
located inside the Bloch ball and its distance to the center of 
the ball is equal to d. 



III. BARYCENTER AS AN ENTANGLEMENT 
MEASURE 

We propose a quantity designed to characterize entan- 
glement for permutation symmetric states. It is based 
on the distance between the barycenter of the Majorana 
points representing a permutation-symmetric state and 
the center of the Bloch ball. We will show that this quan- 
tity vanishes on separable states and it does not increase 
under local operations which preserve the symmetry. 



A. Entanglement measure for permutation 
symmetric states 

For any permutation symmetric state ([3]) with | (fn) = 



cos f | 0} 



sin tH 1) let us define: 



E B = 1 - d 2 , 



(8) 



where d is the distance between the barycenter of the 
Majorana points representing the state and the center of 
the Bloch ball: 



1 

— > (sin 6>i cos $i , sin#isin$i, cos( 



(9) 



where denotes the length of a vector x. 

Let us call Eb the barycentric measure of entangle- 
ment. This quantity can be interpreted as the variance of 
all Majorana points representing the state. Indeed, since 
a vector directed toward any point at the unit sphere has 
the unit length we have E B = D 2 X = (X 2 ) - (X) 2 = 
1 - (X) 2 = l-d 2 . It is easy to see that E B G [0, 1]. 

Note that the quantity Eb can be used to characterize 
the symmetric states of the n-qubit system or the states 
of a single system of size N = n + 1. For a separable 



state |?/>) = \cf))® n (or an 5 , J7(2)-coherent state of a single 
quNit) all stars are located in a single point, so d = 1 
and Eb = 0. Moreover, a generic random state is likely 
to be highly entangled (or highly non-coherent) as the 
Majorana points are distributed almost uniformly at the 
Bloch sphere PHH, which implies d^O and Eb ~ 1 ■ 

A useful entanglement measure should not increase un- 
der local operations. Since the barycentric measure Eb 
is defined only for symmetric states so monotonicity can 
be checked only for local unitary operations which pre- 
serve the permutation symmetry. Consider two n-qubit 
permutation symmetric states \ip), \<j)) connected by a 
local operation. This means that there exist invertible 
unitary operators Ai such that \ip) = Ai ® • • • (g) A n \4>). 
Mathonet et al. proved in [2o| that in fact we may find 
a single invertible operator A for which — A® n \<fi}. 
Geometrically such an operation corresponds to a rota- 
tion of each Majorana points around the same axis by 
the same angle or, equivalently, a rigid rotation of the 
whole Bloch sphere. Obviously thus the radius of the 
barycenter d, and consequently the barycentric measure 
Eb, do not change. This concludes a proof of 

Proposition 1 The barycentric measure Eb is invariant 
under local operations preserving the permutation sym- 
metry. 



B. Examples 

We shall calculate the barycentric measure Eb for cer- 
tain exemplary quantum states and make comparison 
with the geometric measure of entanglement Eq defined 
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FIG. 3. Dependence of entanglement measures for a sym- 
metric 2-qubit state determined by |0i) = |0) and \<j>2) = 
cos ||0) + sin 1 11). Thick curve denotes the barycentric mea- 
sure Eb, while the geometric measure Eg is represented by 
the dashed curve. 
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FIG. 5. Entanglement measures of symmetric 3-qubit state 
of the form such that |0i) = 0, \4>2) = cos ||0) — sin and 
\<p 3 ) = cos f |0) + sin Thick line - E B , dashed line - E G - 
Result for E G from 




ID 



FIG. 4. (color online) Majorana representation of a symmetric 

:os||0) -sin§|l) and \cf> 3 } = 
MP are shown as red, green 



3-qubit state \<f>i) = |0), \(j>2) 
cos||0) +sin||l) with 6 = | 
and blue balls. The black cube denotes the barycenter of MP. 



as a function of the distance from the analyzed state <p 
and the set of separable states T-Lsep- Such a quantity, 
first proposed by Brody and Hughston Q and later used 
in P, can be defined by 



E G {\4>)) = min log 2 



mm 



(10) 



Let us start with the two-qubit case. Without loss of 
generality we may set |0i) = |0) and \(j)2) — cos||0) + 
sin 1 1 1). The special case 9 — |7r is shown in Fig. [2] 
Direct calculation gives for this family of states Eg = 
1 — cos 2 | . Fig. [3] shows how Eb and Eq change with 
the angle 9. We see that in this case Eb > Eq, but for 
a product state (9 ~ 0) and for the maximally entangled 
state (9 = 7r) both quantities coincide. 



Analysis of 3-qubit states exhibits a more interest- 
ing situation. Consider a one-parameter family of 3- 
qubit states, \<fri) = |0), = cos||0) 



sin 1 1 1; 



\<P 



The case 9 



and 

3/ = cus ^|0) + sin||l). The case 9 = |7r is shown in 
Fig. SI Using simple geometry we obtain that for such 
states the barycentric measure reads Eb = 1 — | 2cos 3 9+1 1 2 . 
Fig. [5] shows how Eb and Eq change with the angle 
9. Note a significant difference between Eb and Eq. 
The barycentric measure Eb reaches its maximum at 
9 = |7r, which corresponds to the GHZ state and we 
have E B (\GHZ 3 )) = E G (\GHZ 3 )) = 1. However for 
9 > |7r Eq increases and reaches its maximal value at 
9 = it, at the state \W 3 ) = ^(|001) + |010) + 1 100)). 
Moreover, it is easy to see that the state \GHZ 3 ) is max- 
imally entangled with respect to the barycentric measure 
Eb- for three points at the sphere their barycenter could 
be in the center of the sphere if and only if they form an 
equilateral triangle. 

As a third example we consider the Dicke states | S n ,k) 
(TJ|. They are characterized by k stars at the south pole 
and the remaining n — k Majorana point occupying the 
north pole of the Bloch sphere. 

The length of the radius of the barycenter is thus easily 
calculated as d = I "~ 2fc I. We would like to make the 

1 n 1 

comparison with Eq in a slightly different way employing 
results of Hayashi et al. [2(| who showed that the product 
state closest to | S n ,k) reads 



hence 



|A) = 



E G (S n ^k) = log 2 





(11) 



(12) 



-2k I 



The projection on z axis of A has the length equal | - 
just like the length of the radius d (see above). For 
an even n the state | SVi,n/2) is maximally entangled 
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FIG. 6. (color online) Dependence of entanglement of Dicke 
states | S n ,k), as a function of k for a) n — 10 and b) n — 11. 
Red squares correspond to the barycentric measure Eb, and 
blue circles denote the geometric measure Eg normalized to 
unity for maximally entangled states - lines are plotted to 
guide the eye. 




FIG. 7. Parameterization of the family of the maximally en- 
tangled 4-qubit states. The dashed lines are the 3 axes of the 
coordinate system. The full lines are the 4 sides of the rectan- 
gle. The upper side is parallel to the X-axis, the bottom side 
is rotated by an angle $. The dotted lines are used to mark 
the angles: the bottom dotted line is parallel to the X-axis 
i.e., at the position of the bottom side of rectangle without 
rotation. 



A. Two and three qubits 

For the 2-qubit problem one can establish the following 
property. 



with respect to the measure Eg among all Dicke states; 
for odd n there are two such states, | S , n ,(n-i)/2) an d 
I "Sn,(n+i)/2)< The same states are extremal with respect 
to the barycentric measure. 

Moreover, for an even n the state | S n n / 2 ) belongs to 
the class of the maximally entangled n-qubit states as its 
barycentric measure Eb is equal to unity. A comparison 
between the measures Eb and Eg is shown in Fig. |6]for 
Dicke states with n = 10 and n = 11. By construction 
the same results hold for the states of the standard \j, m) 
basis in the space of N = n+1 level systems, with j = n/2 
and m = k — n/2. 



IV. FAMILY OF MAXIMALLY ENTANGLED 
STATES 

It is interesting to look for states which are maximally 
entangled with respect to the barycentric measure Eb- 
The answer is simple for a small number of qubits. For 
n = 2 a constellation of two stars on the sphere has its 
barycenter in the center of the sphere if and only if it con- 
sists two antipodal points on the sphere. Therefore the 
class of maximally entangled 2-qubit states with Eb = 1 
coincides with all states equivalent to the Bell state up 
to a local unitary transformation. 



Proposition 2 There is no universal unitary operator 
U(6) which transforms any separable state with both MP 
at 9q into a state with MPs at 0q±9 for a given angle 9. 

Indeed, assume that A is such an operator. Then, in 
particular, it would transform the states 1 00) and |11) 
into the Bell state |* + ) = ^(|01) + 1 10) ) , which is im- 
possible since a unitary operator cannot transform two 
orthogonal states into the same state. 

The 3-qubit case was discussed in the preceding section 
IIIIBI with the conclusion that the class of maximally en- 
tangled 3-qubit states consist of the GHZ state and the 
states unitarily equivalent. The case of 4-qubit states re- 
duces to a geometric problem of placing 4 points on the 
sphere in such a way that their barycenter is located in 
the center of the sphere. Consider a rectangle inscribed 
into a sphere, belonging to the XZ plane, so that its 
barycenter is the center of the sphere. Without changing 
the position of the barycenter we can vary the lengths 
of its sides by changing the angle 9 between one side of 
the rectangle and its diagonal (see Fig. [7|). We can also 
rotate one of its sides (the bottom side in Fig. [7]) varying 
the angle $ e [0, tt]. 



B. Four qubits 

Assume now that these four points represent a certain 
4-qubit permutation symmetric state \ip rec }- It means 
that \iprec) is the symmetrization ([3]) of the tensor prod- 
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FIG. 8. (color online) Geometric measure of entanglement 
for the family (|14[) of four-qubit states for which Eb = 1 
parameterized by angles (6, <E>) is not smaller than 1. Cylin- 
der represents the state \Wa), a ball denotes the tetrahedron 
state for which Eg is maximall, while two cubes represent two 
equivalent \GHZa) states (empty cube - \GHZ4) with stars 
situated on the equatorial plane, full cube - rotated \GHZ4), 
situated on a plane perpendicular to the equatiorial one). 



uct of 



\<h) = co S -|0) + e l *sin-|l), 
|0 2 )=co S ^|O)-e < *sin^|l), 
\<h) = cos— ^-|0) + sin— ^-|1), 
\<j> 4 ) =cos^y-^|0) -sin-^— 11), 



(13) 



where 6 G [0, 7r/2], $ G [0, it]. Consequently, 
|Vv ec ) = -s^^e-^ sin 2 ^|0000) + 6e i<& sin 2 6»|1111) + 
[4i cos 9 sin $ - 2(cos 2 + 1) cos $] J2 |0011)) , (14) 



where K is the normalization factor. The state |VVec) de- 
fined above and all locally equivalent states form a class 
of 4-qubit states maximally entangled with respect to the 
barycentric measure Eb- This class contains \Wi) state 
(which we identify with IS42))) obtained for 9 = $ = 0, 
such that all four Majorana points belong to a single line: 
two stars are localized at the north pole and two other at 
the south pole. The choice 9 — ir/2 and $ = ? produces 
\GHZ A ) = ^[|0000) + 1 1111>] represented by stars in 
four corners of a square at the equator. Setting 9 = 7r/4 
and $ = we obtain a state locally equivalent to \GHZ4) 
with stars situated on a plane perpendicular to the equa- 
tiorial one. The tetrahedron state IV'tetr) corresponds to 
8 = arccos(l/\/3) and $ = § . 

To compare both measures of entanglement, following 
[r?) | we computed numerically the geometric measure of 
entanglement Eq for the entire family (fl~4"]) of extremal 



states, for which Eb — 1 (see Fig. [8]). All states from this 
class are characterized by Eq > 1, and the minimum is 
attained for a \GHZ4) state, for which Eg = 1. For a 
\W 4 ) state one has E G {\W 4 )) = log 2 (8/3) « 1.415 0, 
while the maximum is obtained for the tetrahedron state 
and (IV'tetr)) = l°g2 3 ~ 1-585 |17|. This observation 
suggests that the both measures are well correlated, so 
the barycentric measure of entanglement, simple to eval- 
uate, could be used to characterize the state in parallel 
to the geometric measure. 

Note that the states called 'Queen of Quantum', the 
most distant from the set of product states with respect 
to the Hilbert-Schmidt 24j or Bures 2l| distance, are for 
small values of N described by symmetric constellations 
of the Majorana points, so they are maximally entangled 
with respect to the barycentric measure Eb- 



C. Composition of states 

Consider two pure states of a simple, iV-level system, 
and ^2) in Wn- A composition of these states, 
|Vi) © 1^2), defined in g|, 

forms a pure state belonging 
to a 2N — 1 symmetric subspace of the tensored Hilbert 
space, %n ® Hn- Both initial states are described by 
n = N — 1 stars, and their composition, represented by 
2n stars, is defined by adding the stars together. In the 
similar way one can compose the states from two Hilbert 
spaces of different dimensions. More formally, one can 
define the composition by its Husimi distribution Q(a), 
which is given by the product of two Husimi distributions 
representing individual states, 



(15) 



The same idea of composition of two pure states of 
a simple system described by the stellar representation 
can be now used for the symmetric states of several 
qubits. Consider a permutation symmetric n-qubit state 
\tp) G V.f n and another m-qubit state \<j>) G Hf m . Their 



composition, written \tp) \<p) € 



n+m 



represents a 



state of n + m qubits, and is defined by the entire sum 
of n + to Majorana points on the sphere. 

Observe that this notion allows us to write the sym- 
metric Bell state, |*+) = ^(|01) + 1 10>) = |1) |0). 

Introducing two, one-qubit superposition states, |±) = 
^(|0) ± |1)) we see that another Bell state |$_) = 

^(|00) — 1 11)) is equivalent to the composition 



>0 - 



On the other hand, any symmetric separable state of n— 
qubit system can be represented by the composition per- 
formed n—1 times, \<f>, . . . , t/>) = |0) — © |0), as in this 
case all the stars do coincide. 

Take any two figures with the same barycenter. Su- 
perposing them one obtains another figure with the same 
barycenter. This fact implies a useful 

Proposition 3 . Consider two permutation symmetric 
states of any number of qubits, maximally entangled with 
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respect to the barycentric measure, Eb{4>) = Eb(iP) — 
1. Then their composition is also maximally entangled, 
E B (4>Qip) = l. 

To watch this proposition in action consider the com- 
position of two maximally entangled Bell states, | V E'+) 
and |$_). Their composition, |* + ) |$_) = |0) |1) 
|+) |— ) represents the state formed by a square belong- 
ing to the plane containing the meridian of the sphere. 
It is then locally equivalent to the state \GHZ4), repre- 
sented by a square inscribed into the equator, also maxi- 
mally entangled with respect to the barycentric measure. 

The notion of a composition allows us two write down 
a generic random state of an n qubits system 

|Vw) := \fa)®\fa)®---®\<t>n), (16) 

where fa, i = 1, . . . n, denotes one-qubit state, generated 
by a vector taken randomly with respect to the uniform 
measure on the sphere. For large n such a generic state 
is highly entangled with respect to Eb , as the barycenter 
of the points will be located close to the center of the 
sphere. 

Making use of proposition [3] one may also design a 
random state, for which the barycentric measure is equal 
to unity. It is sufficient to compose several maximally 
entangled Bell states represented by a random collection 
of pairs of antipodal points, 

IV4„d) := I0i}0|<M©---0K/2}0K/2>, (17) 

The number n of qubits is assumed to be even and the 
directions labeled by fa and fa are antipodal, so the state 
\fa) \fa) is maximally entangled and locally equivalent 
to the Bell state. Their composition produces the state 
IVrand)' which is characterized by the maximal possible 
value of the barycentric measure, Eb = 1- 



V. DYNAMICS OF MAJORANA POINTS 
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FIG. 9. Dependence of the velocity V = 8/36 of the star 
traveling across the sphere defined as a function of the phase ft 
for the state \fa) = cos^|00) + ^(|01) + |10». The minimal 
V is achieved for a Bell state, /3 = it/2, while the maximal 
velocity is obtained in vicinity of product states, /3 — ► and 

/3 ->• 7T. 




FIG. 10. (color online) Trajectories of Majorana points for 
the state -J=(cos£ | 000) + sin/3 | GHZ}). Regular green, 
dashed blue and thick red lines represent trajectories of 3 MPs 
for /3 £ [0, ir], black dashed line represents the equator. Three 
stars travel from the north pole to the south pole and at the 
equator, for /3 = 7r/2, they form a \GHZ) state. Arrows at 
the red trajectory show the movement direction of the corre- 
sponding Majorana point. 



A. Interpolation between states 

In this section we analyze the parametric dynamics of 
the Majorana points when a state changes with some 
parameter 0. A related problem was analyzed by Prosen 
[27J who analyzed the parametric dynamics of stars on 
the complex plane representing a quantum state under a 
change of paremeters of the system. We will allow the 
pure state to evolve in time, so the parameter /3 can be 
interpreted as time. 

As a first example consider the state \ipi) = cos /3|00) + 
^(|01) + 1 10) ) with P e [0, vr]. It can be written in the 
form ([3]) with the help of the one-qubit states, 



fa) = I 0) 

fa) 



cos/3 I 0) 



K(J3) 



sin/3 I 1), (18) 



where K(f3) 



4 

2-cos^ ■ 



A unitary operator U which 



generates such a dynamics is not uniquely determined. 
If we assume that the operator U transforming |00) into 
\ipi) has the form U = 1 M\ + M2 1, we arrive at two 
possibilities, 



U\ — cos j3 1 er 2 — i —7=- o- y I H 1 a x , 

V2 



U2 — cos P a z I — i S ^ n _ I ' 17 



sin (3 
sin/3 



V2 V V2 



a x I. (19) 



The operators U\ and U2 are similar: the only differ- 
ence is the order of the operators in the tensor product. 
Operating on a permutation symmetric state we can act 
either on the first or on the second subsystem obtaining 
the same outcome. 

Eq. (TT5|) shows that the dynamics of the Majorana 
points for such a family is relatively simple. A single star 
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representing 1 0i) stays at the north pole all the time, 
while the second star makes a loop around the Bloch 
sphere. In particular for (3 = or n we have both points 
at north pole and for j3 = ir/2 we have the Bell state 
^j(|01) + 1 10) ) with one point at the north and the sec- 
ond at the south pole. As shown in Fig. [HI the velocity of 
this second Majorana point is highly nonlinear In the case 
considered the velocity is anticorrelated with the amount 
of entanglement. For the maximally entangled Bell state 
which emerges at /3 — it/ 2 the velocity reaches its min- 
imum while the maximum is attained for the separable 
state 1 00). 

In general, the dynamics of the stars defining the stellar 
representation of a state is not as simple. For instance, 



for a family 



(cos/3 | 000) +sin^ | GHZ)), with 



K(f3) a suitable normalization factor, one obtains the 
motion of 3 points shown in Fig. 1101 



B. Permutation invariant Hamiltonians 

Consider a more general dynamics defined by a one- 
parameter subgroup of the unitary group, 

U = e~ i0H , (20) 

where H is a Hermitian Hamiltonian. Of interest are only 
transformations preserving the permutation symmetry of 
states. We are going to study an exemplary family of 
two-qubit Hamiltonians, parametrized by indices i,j = 
0, . . . 3 with i 7^ j, 

1 



Hi 



- {<Ji ® a j 



where Oi represents the Pauli matrix and ctq = I 
ai<jj = ieijk&k + My we obtain, 

U = -i(cos/3(I<g)]I + cr fe <g)er fe ) + 

i sin/3(cr l ® <jj + a ® Oi) + (I <g> I — at ® Cfc" 

Clearly U preserves the permutation symmetry as it 
acts on both particles in the same way. This is also true 
for any linear combination of generators Hij of the form 
(|2l]>. Taking thus 



(21) 
Using 

(22) 



H = 



H 



23 



H, 



1 



V2 



2 v / 2 



1 (cr + cr 3 ) + (cro + er 3 ) ® a 2 



we obtain: 



U = 



I cos (3 

sin 3 

~7T 

sin d 

V 



sin j3 

i(l + cos/3) 
i(-l + cos^) i(l + 




V2 \ 

i(-l + cos/3) 

i(l + cos/3) 

1/ 



(23) 



which transforms the state |00) into \4>i) = cos/3|00) + 
^^(|01) + 1 10)) as in the example considered previously. 
Note that the operator U which generates the same dy- 
namics of stars differs from the operators U± and U2 de- 
fined in (HI 



C. Unitary dynamics of symmetric states of (TV — 1) 
qubits and corresponding dynamics of a single quNit 



As discussed in previous sections, any symmetric pure 
state of n-qubit system can be associated with a state 
of a simple system of size N = n + 1. Thus any unitary 
dynamics of the composite system, which preserves the 
permutations symmetry, induces a certain dynamics in 
Hn- To make such a link explicit consider a permutation 
symmetric state if) € 7if n which can be represented in 
the standard basis, {|0)® fc ® \l)^ n ' k) } <k<n- Define a 
a unitary transition matrix T that transforms this basis 
into an orthonormal basis with first n + 1 permutation 
symmetric vectors and remaining 2™ — (n + 1) vectors, 
which span a basis in the orthogonal subspace. Denote 
|V/) = T~ l \il>) = (g) where \4>) £ H n +i, |0> € U 2 n- n -l 
and take any unitary operator U, acting on T-Lf n , which 
preserves the permutations! symmetry. Wc have then: 



T _1 |-0) = T^UTT^ 1 ] 

W) = u'W), 



(24) 



where U' = T~ X XJT and \tp') = T^ 1 ^). As U preserves 
the permutation symmetry condition Eq. (|25[) implies the 
following block structure of the matrix U': 



U' 



V 








w 



(25) 



Here V and W denote unitary matrices of size n + 1 and 
2 n — (n + 1), respectively. The matrix V acts on the 
(n + l)-dimensional subspace of permutation symmetric 
vectors. These vectors are in one to one correspondence 
with vectors spanning an orthonormal basis in Tin+i ~ 
see e.g. (28|. In other words the matrix U representing 
a permutation symmetry preserving dynamics has to be 
reducible. The unitary block V of the rotated matrix, 
U' = T~ 1 UT : defines thus a unitary dynamics of the one 
quNit system, associated with the n-qubit dynamics U, 
which preserves the permutational symmetry. 



D. A three— qubit example 



To illustrate the above reasoning let us ana- 
lyze a 3-qubit example. We use the natural 
basis, (|000), |001), |010), |011), |100), 1 101) , |110), | 111)) 
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and the unitary transition matrix T reads: 



(I 


























1 








1 


1 

76 











1 

75 








1 


1 

\/8 














1 











1 

75 


1 

75 





1 











V 3 














1 











1 

V2 


1 








1 

7a 














Pi 

V 3 


U 








1 















(26) 



/ 

Using formula (1201) we construct permutation symmetry 

'lo\ 

v o o y 

where ^ stands for the sum over all permutations. The 

7T 

matrix U obtained in this way does not have any special 
structure. However, the transformation T brings it to the 
matrix U' = T~ 1 UT which enjoys the block structure as 
in (US]), 



preserve operator U taking H = ^2 a\ ® 03 



V = 



i(4<3) 



■ JiV3»in(4/3) 2v / 3cos(,3) 2 si,,(,3) 2 
cos(43) ^isin(4/S) 
ii.ir»(40) 1(3 + cos(4/3)) 



The matrix V determines thus the corresponding dy- 
namics of a simple system of size n + 1 = 4, while the 
remaining part, 




FIG. 11. (color online) Trajectories of MP for the state 
\ipn) = cos/3 I 00) — sin/3 | 11). Red and green line repre- 
sent trajectories of two points for /3 £ [0, n], black dashed line 
represents the equator. Both stars travel from the north to 
the south pole, where they meet at /3 = n/2. Later on they 
return back to the north pole. 



W 



co.03) 




is irrelevant for the corresponding dynamics in 



FIG. 12. (color online) Velocity 9 for the state l^) = cos/3 | 
00) — sin p I 11). Full line is 9 in the function of P, dashed 
line is 9{j3). 



E. Dynamics for two qubit case 



From (f20|) we obtain a differential equations for the 
dynamics 



id/3 \ip)=H \ V). 



(27) 



which can be also written in terms of the angles 9i and 
the phases $i, i = 1,2 characterizing the positions of 
the Majorana points. In the general case such equations 
are not easy to solve. Therefore we will restrict here our 
attention to the special case of n = 2 qubits. 

Let us take H = —o\ ® 02 — 02 ® o\ . Using Eq. ([20|) 
we obtain the operator U which generates dynamics of 
the MPs of state |t/> 2 > = cos ,3 | 00) - sin/3 | 11). The 
corresponding dynamics of the stars on the sphere with 
respect to the time j3 is shown in Fig. [TT1 

For simplicity we restrict our analysis to j3 € [0,7r/2], 
where B\ = Q%, $1 = $2 + tt and the phases $1 and $2 



do not change. Using Eq. (|2~T1) we obtain a differential 
equation for the polar angle = 6\ = 8 2 , 



3 + cos 36* 
2 sin 9 



(28) 



As shown in Fig. [T2] the velocity diverges when 9 ap- 
proaches to zero or 7r and both stars are located in a 
single point. The velocity 9 reaches its minimum for a 
Bell state obtained for 9 = n/2 

It is instructive to compare the velocity of stars during 
the unitary dynamics of a state parametrized by the the 
phase ft with changes of the barycentric measure Eb - 
see Fig. [T31 Observe an anticorrelation between entangle- 
ment measure and 9 - the maximal velocity occurs for the 
minimal entanglement and vice versa. This observation 
implies that product states tend to gain entanglement 
during a relatively small interaction time. 
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0.0 0.5 1.0 1.5 



9 

FIG. 13. Dependence of the barycentric entanglement mea- 
sure Eb on the phase ft for the state \^2) = cos ft | 00) —sin ft \ 
11). 



VI. CONCLUSIONS 

In this work we presented a homogeneous approach 
to study the structure of pure quantum states describing 
two physical problems: a simple system consisting of n+1 
levels and the class of states of an n-qubit system, sym- 
metric with respect to permutations of all subsystems. 
Making use of the Majorana-Penrose representation one 
can find a direct link between these two cases as any con- 
stellation of n stars on the sphere determines a quantum 
state in both setups. In particular, product states of the 
multi-qubit system, for which all Majorana points coa- 
lesce in a single point, correspond to spin coherent states 
of the simple system. To find a direct relation between 
two problems one can identify the orthogonal basis of the 
eigenstates of the angular momentum operator J z acting 
on rl n +i with the set of (n + 1) Dicke states which span 
the complete basis in the subspace of symmetric states 
of the composite system. 

Physical properties of a given pure state can be thus 
related to the distribution of the corresponding collec- 
tion of n points on the sphere. The variance of this dis- 
tribution, related to the radius of the barycenter inside 
the ball, can be thus used to characterize the degree of 
non spin coherence of the states of a simple system or 
the degree of entanglement for the composite systems. 
The proposed barycentric measure of quantum entangle- 



ment achieves its maximum for these states, for which 
the barycenter of the corresponding Majorana points is 
located at the center of the ball. This class of states in- 
cludes the Bell state of a two-qubit system, the GHZ state 
of a three qubit system and several states distinguished 
by being most distant from the set of separable states 
and called 'Queen of Quantum' In the case of four- 
qubit states we have explicitly described a two-parameter 
class of extremal symmetric states for which the barycen- 
tric measure achieves the maximal value, Eb = 1- All 
these states are also highly entangled with respect to the 
geometric measure Eq, which for them belongs to the 
interval [1, log 2 3]. 

It is convenient to define the composition of two per- 
mutation symmetric states of an arbitrary number n and 
m qubits, each of them described by n and m Majorana 
points, respectively. The composition is characterized by 
the collection of n + m stars on the sphere and represents 
a symmetric state of m + n qubits. Making use of these 
notion we show that the state \GHZ4} can be interpreted 
as a composition of two two-qubit Bell states and con- 
struct a family of maximally entangled multiqubit pure 
states with Eb = 1. 

Any unitary dynamics acting on the n-qubit system is 
described by a matrix U of order 2". Under assumption 
that the dynamics does not break the permutation sym- 
metry, the matrix U is reducible and can be written as 
a direct sum of two unitary matrices, U = V ffi W. The 
matrix V or order n + 1 describes the unitary dynamics in 
the subspace of symmetric states of the composed system 
or the corresponding unitary dynamics in the space of all 
states of the (n + 1) -level system. A unitary dynamics of 
a quantum pure state leads to a non-linear dynamics of 
the corresponding stars of its Majorana representation. 
In a simple model evolution investigated for a two qubit 
system the velocity of stars is small if they are far apart, 
what corresponds to the highly entangled states, and it 
increases as the stars get together and the state is close 
to be separable. 
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number N N202 261938 of Polish Ministry of Science and 
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